Earth impact threat of newly discovered 
asteroids: an a posteriori analysis 



Germano D'Abramo 

INAF - Istituto di Astrofisica e Planetologia Spaziali, 
Via Fosso del Cavaliere 100, 00133, Roma, Italy 
E-mail: Germano . Dabramo@inaf . it 

a 

h— 5 Abstract 

In the past decade both scientists and laymen have probably heard at 
s -! least once through the newspapers, TV and Internet that a new aster- 

oid has been discovered with non-zero (sometimes "high") probability of 
' pH ' collision with the Earth in the near future. Since early 2000's, such prob- 

er | abilities are routinely calculated by two impact monitoring systems (one 

I in the US, the other in Europe) on preliminary orbits of newly discovered 

Near Earth Asteroids (NEAs) , and are regularly updated whenever addi- 
tional new astrometric observations for each threatening object become 
Q | available. A typical pattern is that as the orbit becomes more precisely 

CZ3 determined, impact probability often increases initially, but then turns 

around and decreases until it falls to zero, or some very low number. 
O In the present study we define a probability measure which provides a 

' ^5 simple tool to evaluate from the very beginning the chance that the im- 

pact probability calculated by monitoring systems for a threatening object 
will reach unity at the end of the overall process of orbit refinement. We 
stress that this chance is independent of the specific (fluctuating) value of 
the probability calculated by the monitoring systems. In a precise sense, 
which should be clear throughout the paper, the concrete impact proba- 
bility of a newly discovered asteroid is not that given by the monitoring 
systems. 

(N 
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O 1 Introduction 

Starting in the 1980s, news and announcements about newly discovered Near 
Earth Asteroids (NEAs, see Fig. [I] for the astronomical classification of such 
planetary objects) with unusually high Earth-collision chances in the near future 
are appearing periodically on the press, TV and recently on the Web. 

Admittedly, over time there is dwindling feedback of public opinion, unless 
the announcement concerns a new, relatively big asteroid (hundreds of meters in 
size) with an alarming high probability of hitting the Earth in the next decades. 

To date, impact monitoring and analysis are carried out exclusively by two 
research groups in the world, one at the University of Pisa (the system is called 
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Figure 1: Near Earth Asteroids, or NEAs, are all those asteroids with perihe- 
lion distance q less than 1.3 A.U. (1 A.U. is the mean Earth-Sun distance and 
it is nearly equal to 1.5 x 10 s km). The perihelion is the point in the orbit of a 
planetary object where it is nearest to the Sun, as opposed to aphelion where 
it is farthest. NEAs are divided in three sub-classes: Apollos, Amors, Atens. 
Apollos have q < 0.983 A.U. and semi-major axis a > 1A.U., Amors have 
1.017 A.U. < q < 1.3 A.U. and Atens have aphelion distance Q > 0.983 A.U. and 
semi-major axis a < 1 A.U. Note that 0.983 A.U. and 1.017 A.U. are the perihe- 
lion and aphelion distances of the Earth, respectively. The above picture shows 
a projection of NEA orbits onto the ecliptic plane (the Earth's orbit plane). 
That could give the wrong impression that actually all Apollo and Aten orbits 
intersect the orbit of the Earth. Obviously, this is not true since orbits are also 
inclined with respect to the ecliptic plane, and rotated around the focus (the 
Sun) 



CLOMONQ and the other at NASA-JPL (where it is called SENTHVQ. They 
began routine operations in the early 2000's. On a daily basis, these two groups 
collect all astrometric data of newly discovered NEAs; actually, they collect 
data of all NEAs, new and known, observed the previous night by the observing 
surveys. 

In this study we do not take into account the impact threat due to cometary 
bodies. Although comets (especially long-period comets with orbital periods in 
excess of 200 years) are believed to present less than of the order of 1 percent 
of the total Near Earth Objects (NEO) impact threat |Stokes et al. (2003)| , 



their contribution to impact flux and threat to the Earth cannot be neglected 
a priori, as demonstrated by the extraordinary impact of Comet Shoemaker- 
Levy 9 on Jupiter in July 1994. Nevertheless, impact monitoring systems do 
not currently handle cometary bodies for a variety of reasons, among which: the 
greater orbit uncertainty of comets due to the erratic non-gravitational forces 
which affect them and their lower intrinsic risk due to higher orbital periods 
(which are often higher than the 100 year time-limit in the monitoring systems' 
current predictions). 

Astrometric data, or astrometric observations, are the optical/radar mea- 
surements of the asteroid sky position with respect to the so-called fixed stars 
(the Earth, in the case of radar observations). These data come from all the 
observatories around the world and are publicly available at the Minor Planet 



: http : //newton . dm.unipi . it/neodys/index.php?pc=4. 1 



http : //neo . jpl.nasa.gov/risk/ 
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Center (MPCf] at Harvard, Massachusetts. SENTRY and CL0M0N2 then 
daily apply their algorithmic procedures in order to identify possible threat- 
ening asteroids among the new ones and to refine the orbit of those already 
categorized as threatening. 

Probably, the first well-publicized scare concerned asteroid 1997 XFn, for 
which a possible impact with the Earth was predicted in 2028 by Brian Marsden 
at MPC. This object and its announced non-zero impact probability in the near 
future rose a sort of unprecedented turmoil within the astronomical community 
and, to a lesser extent, also within the public opinion. Instead, one of the 
first case in the history of the impact monitoring era (that of CLOMON2 and 
SENTRY) was that of asteroid 1999 AN 10 : this object played an important role 
in developing the concept of Virtual Impactors described later in this paper (for 
more information, see Milani et al. (1999) [Milani et al. (2003)] and the News 



Section of the NASA Ames Research CenteiFl for a detailed historical account 



of these cases, see [Chapman (1999)] and [Marsden (20"04)] ). 



For the record, one of the latest asteroid discovered with near future non- 
zero impact probabilities is the ~ 180 meters sized 2011 UE305, while probably 
the most famous one is asteroid (99942) Apophi^ (~ 300 meters in size): in 
December, 2004 monitoring systems calculated an initial impact chance for the 
year 2029 as high as ~ 1/38 [Sansaturio and Arratia (2008)| . Then, after more 
observational data were gathered, the impact in 2029 was ruled out. Currently, 
Apophis has a chance nearly equal to 1/230.000 of hitting the Earth in 2036. 
For a daily updated list of all potential future impactors, see CLOMON2 and 
SENTRY Web pages. 

The discovery of new NEAs with non-zero impact probabilities always rises 
the attention of the astronomical community, at least of those who work in 
that field, and sometimes, as said before, the attention/interest of the general 
public: this typically happens when the calculated impact probabilities for a 
specific date in the future result to be unusually high according to the standards 
of the astronomical community (see the International Astronomical Union - 
IAU- rules^J and when the researchers stress that the case should deserve more 
attention. 

To be more precise, unusually high means that the impact probability cal- 
culated by monitoring systems for a particular new asteroid is higher than the 
yearly statistical impact probability on the Earth coming from an unspecified 
asteroid of equal size or larger belonging to the whole NEA population^] multi- 
plied by the time interval (in years) which separates the present time from the 
future impact date of the asteroid under analysis (see Chesley et al. (2002)] ). 



In short, a newly discovered asteroid with non-zero chance of hitting the Earth 
on a specific date in the future deserves attention if its calculated impact proba- 
bility is greater than the statistical probability that an unknown asteroid of the 
same size or larger is likely to hit the Earth over the years until that date. Ba- 



3 http : //www. minorplanetcenter . org/iau/mpc .html The MPC is responsible at inter- 
national level for the designation of minor bodies in the Solar System and for the efficient 
collection, computation, checking and dissemination of astrometric observations and orbits 
for minor p lanets and com ets. 

^http : //impact . arc .nasa. gov/news_detail . cf m?ID=48 

5 It is named after an ancient Egyptian evil god. 

e http : //web. mi t . edu/rpb/wgneo/Te chComm.html 

' Due to the natural fact that our planet moves in an environment swept by the NEA 
population. 
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Figure 2: Cumulative population of NEAs vs. absolute magnitude, with 
equivalent scale for diameter in km, and expected impact interval in years on 
the right (this plot is a modified and updated version of a similar figure in 
Morrison et al. (2003)] ) . The various points represent estimates of NE A popu- 
lation, and thus estimates of the Earth impact flux, by different methods (whose 
description is beyond the scope of the present paper). The straight line is a 
simple power law that approximates the data: with respect to the right hand 
vertical scale, this line represents j-j^jyj, where pi{> D) is the equation nil) 
of the text. The two downward arrows indicate the Tunguska and K-T (end- 
Cretaceous) impacts. The dot-dashed line is the cumulative NEA population 
discovered to 06/30/2006 (courtesy Alan W. Harris) 



sically, this is the philosophy of the hazard scale most used among astronomers, 
the Palermo Scale |Chesley et al. (200 2) , where an impact prediction (time, 
size and probability) is normalized to the background expected impact probabil- 
ity from now to the time of the prediction. Expressed as a logarithm, Palermo 
Scale 0.0 means the predicted event is essentially "as expected" in that time in- 
terval; positive values mean the predicted event is "extraordinary" and deserves 
attention; and negative values indicate the predicted event is only a small ad- 
dition to the expected impact flux. 

The background impact probability is extrapolated from a statistical analysis 
of NEA close encounters with the Earth and from the lunar cratering record 
Morrison et al. (2003) ). It be fairly approximated in the following 



see 



exponential form: 
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M>^) = 10£-"yr~\ (1) 

where D is the diameter of the asteroid expressed in meters. Rigorously speak- 
ing, it gives the average number of NEAs larger than a given size that hit the 
Earth per year. If this number results to be less than unity, as happens with 
bigger asteroids, then it properly expresses a mathematical probability; other- 
wise, if this number results to be greater than or equal to unity, as happens 
with (sub) meter-sized NEAs that fall on the Earth every year more than once, 
then it more properly expresses a kind of frequency of impacts. Accordingly, the 
average number of NEAs that hit the Earth per year with diameters between D 
and D + AD is given by: 

Pi(D + AD; D) = 10(L>~ 2 - 3 - (D + AD)" 2 3 ) yr" 1 . (2) 

The reciprocal of the background impact probability ( p .r>jj) ) is the impact 
interval time and it gives the mean time between two consecutive impacts of 
asteroids larger than a given size (Fig. [2]). 

As is clear from Fig. [2] the constant power law ([I]) is only a bona fide ap- 
proximation of various data coming from different sources: the points plotted 
in Fig. [2] represent the estimate of the NEA population (and thus, the estimate 
of the Earth impact flux through the statistical analysis of close encounters 
with the Earth of the estimated NEA population) in different size ranges and 
by different methods. It must be said that every estimate has its own intrinsic 
uncertainty which is often very high and sometimes it cannot even be properly 
evaluated. By the way, the central part of the plot comes from the lunar crater- 
ing record and from a mathematical debiasing of the observational survey data, 
while the flux of meter-sized objects (upper left of the Figure) is obtained from 
detections of small objects entering Earth's atmosphere Silber et al. (2009)] . 



It is clear that the background impact probability is anything but well- 
constrained; however, this is not a concern for us, since the idea at the basis of 
our argument is somewhat independent of the background impact probability 
exact form. Moreover, eq. is the same commonly used within the planetology 
community, e.g. in the definition of the Palermo Scale. 

Returning to the impact monitoring procedures, if and when additional new 
astrometric observations become available, the asteroid orbits are refined and 
usually the impact possibilities are definitively ruled out (more on this in Sec- 
tion 2), as happened with the case of asteroid 1999 ANio- Since that time, the 
interest on the single case typically returns to be a purely academic one. 

In the following Section we qualitatively describe how the current impact 
monitoring systems compute impact probabilities of newly discovered NEAs. 
In the third Section we define the a posteriori conditional probability W: it is 
strictly related to the background impact probability and to the mean annual 
frequency with which impact monitoring systems find impact threats among 
newly discovered asteroids. We argue that probability W is the proper proba- 
bility measure which helps to assess from the very beginning the future impact 
chances of asteroids just inserted in the risk lists of monitoring systems. Of W 
we also give an upper-limit estimate. 
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2 NEA discovery and impact scare 



Soon after a new NEA is discovered by telescopic surveys around the world, 
a preliminary orbit is computed using its positions in the sky over a suitable 
(minimal) interval of time (astrometric observations). Like every physical mea- 
surement, astrometric ones are affected by errors which make the resulting orbit 
uncertain to some variable degree. Sophisticated mathematical and numerical 
tools are now available to orbit computers which allow to propagate such mea- 
surement errors to the six orbital elements which identify the orbit of the aster- 
oid (semi-major axis a, eccentricity e, inclination i, longitude of ascending node 
f2, argument of perihelion cj, time of perihelion passage t). For such reasons, 
the new NEA, soon after its discovery, is not represented by a single point in 
the 6-dimensional orbital elements space; rather, it is represented by an uncer- 
tainty region, a 6-dimensional volume with blurred contours. Obviously, the 
volume of this uncertainty region changes (usually it shrinks) when additional 
new observations become available and the orbit refines. 

Moreover, when the nominal orbhj^Jof the new NEA is geometrically close to 
the orbit of the Earth, and it shares some other peculiar orbital characteristics 
(like encounter timing issues), it can happen that some orbital solutions which 
lead to a future collision of the asteroid with the Earth can not be excluded only 
on the basis of the available astrometric observations. Namely, orbital solutions 
which lead to a collision are inside the uncertainty region and they are fully 
compatible with the available astrometric observations and their errors. 

What is substantially done in these cases by monitoring systems, with var- 
ious techniques whose description is well beyond the scope of this paper (see 
Milani et al. (2000) | |Milani et al. (2003)| |Milani et al. (2005)] ), is to sample 



the uncertainty region with an appropriate number of sample points according to 
a suitable 6-dimensional space distribution, closely related to what is currently 
known about error statistics (a priori uniform or Gaussian; see Carpino et al. (2003)] ). 



Then they evaluate the relative probability that the 'true' orbit of the asteroid 
is one of the collision ones. Henceforth, we will refer to this probability with 
the symbol Vi. The collision orbits are nowadays commonly called Virtual Im- 
pactors (or Vis; for an exhaustive review see Milani et al. (2003)] ). Sometimes 
we will call Vi as VI impact probability. 

If a newly discovered NEA exhibits Vis, then it is suddenly inserted by mon- 
itoring systems into their publicly available risk lists, together with its proba- 
bility V t and Palermo Scale rating. Currently, CLOMON2 and SENTRY find 
hundreds of newly discovered NEAs with VI orbital solutions every year. For a 
statistics of Vis identification between calendar years 2004 and 2009, see Fig. [3] 

Every time additional new astrometric observations become available, the 
quality of the asteroid orbit improves and the estimated impact probability Vi 
is re-computed (i.e. it changes). This usually happens in the weeks, months 
and even years following the discovery date. A typical pattern is that as the 
orbit becomes more precisely determined, impact probability Vj often increases 
initially, but then turns around and decreases until it falls to zero, or some very 
low number. The reason of the initial increasing behavior is rather technical and 
it is essentially tied to the fact that uncertainty region shrinks with additional 
new observations, while usually some VI orbital solutions still remain inside the 



8 The orbit which fits the observations best. 
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Figure 3: NEAs discovered between calendar years 2004 and 2009 exhibiting 
VI orbital solutions (taken from SENTRY risk list and archive). The increasing 
trend with time is mainly due to the fact that discovery surveys have become 
more efficient in discovering NEAs over the years. The efficiency of impact 
monitoring systems in finding Vis is very high and it has been almost always 
the same through the years 

uncertainty region in the elements space. 

Given the past experience of public and professional reactions mentioned 
in Section I and given the current high rate of VI orbital solution discoveries, 
some questions might arise: how are present VI impact probabilities actually 
related to the concrete future impact? How much threat should we reliably 
read in a Vis detection announcement? In the following Section we propose an 
a posteriori conditional reading of VI impact probabilities. The aim of this new 
probability measure is to provide a consistent way to cope with such questions. 

Within the paper we label our probability as ' a posteriori' since it is obtained 
as the ratio between two statistical quantities (relative frequency). On the 
contrary, the VI impact probabilities can be considered a priori, in the sense that 
they are obtained through sophisticated mathematical models and deductive 
reasoning (for instance, the choice of an a priori 6-dimensional space distribution 
to sample the uncertainty region). 

3 Statistical reading of V^: the probability W 

The concrete question about the future impact chance of a newly discovered 
NEA is the following: soon after the discovery of VI orbital solutions of a new 
asteroid, what is the probability that Vj approaches and eventually reaches 
unity (within this paper we will use the compact notation 'Vj — > 1'), after the 
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right amount of additional new astrometric observations has become available? 
Actually, what we are asking is: by only knowing that a newly discovered NEA 
exhibits some VI orbital solutions, what is the probability that Vi will be equal 
to 1 at the end of the whole orbital refinement process? 

Now, let us make the following thought experiment, only functional to the 
presentation of our point. Suppose that the existing discovery surveys are able 
to discover all NEAs which pass close to the Earth down to a size cut-off. This 
is obviously not true since during their close approaches to the Earth, a lot 
of still unknown asteroids remain too dim to be detected by telescopes: they 
are too small in size and/or still 'too distant'. Moreover, some NEAs are not 
found because telescopic observations miss them, namely surveys do not image 
a portion of the night sky when they are there and bright enough to be seen. 
This becomes especially true for asteroids relatively close to the Sun, where 
monitoring is more sporadic or even impossible. 

We also reasonably suppose that every discovered asteroid really impacting 
the Earth in the future will show some Vis, with low Vi soon after the discovery 
and fluctuating with an increasing trend as soon as subsequent astrometric 
observations become available, as usually happens in reality. 

In other words, we are putting ourselves in the somewhat idealized situation 
where every impacting asteroid above a size cut-off is surely discovered and for 
it monitoring systems surely spot some Vis soon after its discovery. 

Thus, we define W, the a posteriori conditional probability of Vi — > 1, which 
could also be interpreted as a kind of 'weight' of the VI impact probability 
calculation, as: 

W{D + AD;D) = lim <D + ^D;D) 
V ' ; t->oo v(D + AD; D) 

n(D + AD; D) 1 Pl {D + AD;D) 

t^o T v(d+ad-d) f Vt (D + AD;DY 1 ' 

where n(D+AD; D) is the number of asteroids with size between D and D+AD 
which actually impact the Earth in the period of time T, with T> 1 year, and 
v(D + AD; D) is the number of asteroids found by monitoring systems among 
all the newly discovered ones to exhibit VI orbital solutions, in the same size 
interval and in the same period of time T. 

The limit lim-j^oo n ( D+ ^ D > D ) j s ^ e definition of the background annual im- 
pact probability, equation (2h. The function /^(D+AD; D) = lim*r->oo V ( D+ ^ D > D ) 
is the annual frequency ofnewly discovered NEAs with sizes between D and 
D + AD found with VI orbital solutions. 

Note that, according to what we have said at the beginning of this Section, 
the number n(D + AD;D) is counted in the number v{D + AD;D), since we 
have assumed that every impacting asteroid is identified soon after its discovery 
as having some VI orbital solutions, thus n(D + AD; D) is always less than or 
equal to v(D + AD;D). 

Let us still focus our attention on eq. |3]) . Within the hypotheses introduced 
above, we imagine to be able to wait for a very long period of years, T, and count 
the number n(D + AD; D) of true asteroid impacts on the Earth and the number 
v(D + AD; D) of newly discovered NEAs found by monitoring systems to have 
VI orbital solutions with size between D and D + AD and within that period 
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of time. Hence, the limit for T — > oo of the ratio between these two numbers 
gives the a posteriori conditional probability that a discovered asteroid with size 
between D and D + AD, for which monitoring systems have spotted some VI 
orbital solutions, is just that which will surely fall on the Earth. In the fourth 
member of eq. (|3| we have simply rewritten W in terms of the background 
annual impact frequency pi and the mean annual frequency fy i . 

Thus, inevitably W(D + AD; D) gives the probability that an asteroid with 
Vis, and thus an asteroid inserted into the risk lists of the monitoring systems, 
has its probability Vi eventually reaching unity. 

Before we said that W can be interpreted as a kind of 'weight' of the VI 
impact probability calculation. Suppose that, thanks to improvements in ob- 
servational techniques (e.g. higher positional precision) and orbits computation, 
the number of newly discovered asteroids identified by monitoring systems as 
potential impactors sensibly decreases in every diameter (or absolute magni- 
tude) range. Accordingly, the probability W will increase (given its definition) 
in every diameter (or absolute magnitude) range. Since the decrease of the 
number of potential impactors among the new asteroids means an increased 
capability in constraining the true potential impactors, the consequent increase 
of W would equivalently mean a increased capability in constraining the true 
potential impactors by the monitoring system. Thus, W could be seen as a 
'weight' in expressing the actual capabilities of the monitoring system. 

It is also easy to note that the numerical value of W(D + AD; D) is naturally 
bounded as follows: 

N(D + AD;D) ~ y ' ' ~ K ' 

where N(D + AD; D) is the actual number of all NEAs with size between D and 
D + AD. The lower bound derives from the fact that the number v (D + AD; D) 
can not be greater that the entire population of existing NEAs in the same size 
range. For the sake of the argument, in eq. Q we have left aside the homogeneity 
of physical dimensions. 

Moreover, we can see that W is not directly related to the specific numerical 
value of Vi, no matter how Vj's specific, fluctuating numerical figure is. Rather, 
it depends upon fy i which, in turn, depends upon some contingent observational 
characteristics. These features are the annual number of NEA discoveries, the 
amount of astrometric observations available at discovery, the magnitude of 
astromctric errors and the conventions in their statistical treatment, as well as 
the observational geometry and orbital characteristics of the newly discovered 
asteroid. 

Although the value of /y 4 depends upon contingent, variable features, it 
is worthwhile to provide a statistical estimate of it. Given the total number 
of NEAs with Vis found at every size between calendar years 2004 and 2009 
(Fig. [3]) an indicative estimate of fv t {D + AD; D) is possible, see Fig. |4j and 
consequently a preliminary estimate of W(D + AD; D) can be obtained as in 
Fig. 

Note that relaxing the optimistic assumptions on the almost perfect NEA 
discovery efficiency and VI monitoring capabilities makes /y 4 , as approximated 
with the aid of Fig. [4j even an underestimate. And consequently, our computa- 
tion of W(D + AD; D) is surely an overestimate. 
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Figure 4: Annual frequency fvi(D + A; D) of VI detections, estimated using 
all NEAs discovered between the calendar years 2004 and 2009 and listed in the 
SENTRY risk list (and archive) . The size range of this group is from ~ 3 m up 
to a maximum of ~ 4 km. The size bins are spaced geometrically (such that 
D + AD = 1.3D) 



Furthermore, /y 4 changes over time, since the discovery completion increases 
and the discovery rate declines with time (the NEA population is finite). Thus, 
it is useful to update /y ; (and hence W) from time to time. But again, we stress 
that the sense and the validity of the definition of W is not affected by (and 
dependent on) such a time dependence. 

Summarizing, the function W(D + AD; D) provides us with a first simple 
tool to evaluate from the very beginning the chance that Vi, i.e. the impact 
probability calculated by the monitoring systems for newly discovered NEAs 
with VI orbital solutions, eventually reaches unity. As a matter of fact, Vi is 
a stochastic variable since nobody knows how Vi will evolve with additional 
observations, and it is perfectly legitimate and valuable to define a probability 
(W) of a (stochastic) probability value (Vi). 

Consider now the application of probability measure W to two well known 
cases: that of Apophis and that of asteroid 2008 TC3, which actually impacted 
the Earth the day after its discovery. These applications suggest the reliability 
of probability W in providing an early direct glimpse of the most likely fate of 
probability Vi. 

Soon after Apophis was discovered in December 2004, impact monitoring 
systems identified multiple Vis orbital solutions and obtained an initial scary 
impact probability of ~ 1/38 for the year 2029 Sansaturio and Arratia (2008)] . 

But, what is the probability W for Apophis? Apophis is an asteroid with an 
estimated diameter of nearly 300 m. According to Figure [5j W(> 300 m) is less 
than 10 -6 , namely more than four orders of magnitude lower than that initially 
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Figure 5: The probability W(D + AD; D) as a function of asteroid diameter. It 
is thought that the whole population of NEAs with diameter larger than nearly 
4 km has been probably already discovered by astronomers and the orbits are 
also known to be safe for the Earth. This is the reason why no VI has been 
found so far for such NEAs. The size bins are spaced geometrically (such that 
D + AD = 1.3D) 



reported by monitoring systems and close to their current estimate, obtained 
after some orbit refinement. 

Nowadays, it is almost certain that the probability Vi for Apophis will go 
to zero with future astrometric observations, but it was not so clear at the 
beginning of the impact monitoring process, with an initial Vi as high as ~ 1/38. 
As a matter of fact, an early use of probability W would have given a direct 
glimpse of what would have been the most likely fate of Vi for Apophis. 

The same would have happened for 2008 TC3. This asteroid was discovered 
on October 6, 2008 and impacted the Earth about 20 hours later. The impact 
was predicted by monitoring systems to have a probability of ~ 100%. Con- 
sistently, the probability W for 2008 TC3 is very close to unity, the size of the 
asteroid being estimated to be between 2 and 4 m (see Fig [5]). 

4 Conclusions 

In this paper we have defined the a posteriori conditional probability W and dis- 
cussed its properties. W is strictly related to the background impact probability 
Pi and to the mean annual frequency with which the impact monitoring systems 
find Virtual Impactor orbital solutions among newly discovered asteroids. 

This probability measure provides us with a simple tool to evaluate from the 
very beginning the chance that the probability Vi, namely the impact probability 
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calculated for a newly discovered object by monitoring systems, will reach unity 
at the end of the whole orbit refining process. In short, W is a probability of a 
stochastic variable (the probability Vj). Hence, the proper impact probability 
of a newly discovered asteroid is not Vj, which actually fluctuates with new 
additional observations, but the probability that V$ — > 1, i. e. W. 

The application of probability measure W to the cases of Apophis and 
2008 TC3 suggests its reliability in providing an early direct glimpse of the 
most likely fate of probability Vj. 
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